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Abstract: We demonstrate the phase fluctuation introduced by oscillation 
of scattering centers in the focal volume of an ultrasound transducer in 
an optical tomography experiment has a nonzero mean. The conditions 
to be met for the above are: (i) the frequency of the ultrasound should be 
in the vicinity of the most dominant natural frequency of vibration of the 
ultrasound focal volume, (ii) the corresponding acoustic wavelength should 
be much larger than £*, a modified transport mean-free-path applicable for 
phase decorrelation and (iii) the focal volume of the ultrasound transducer 
should not be larger than 4 — 5 times (^) 3 . We demonstrate through sim- 
ulations that as the ratio of the ultrasound focal volume to (^*) 3 increases, 
the average of the phase fluctuation decreases and becomes zero when the 
focal volume becomes greater than around 4(£*) 3 ; and through simulations 
and experiments that as the acoustic frequency increases from 100 Hz to 
1 MHz, the average phase decreases to zero. Through experiments done in 
chicken breast we show that the average phase increases from around 110° 
to 130° when the background medium is changed from water to glycerol, 
indicating that the average of the phase fluctuation can be used to sense 
changes in refractive index deep within tissue. 

© 2012 Optical Society of America 
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1. Introduction 

Ultrasound modulated optical tomography (UMOT) combines the advantages of an optical 
property recovery in a scattering medium with the high resolution available with focused ultra- 
sound (US) [1,2]. Focused US creates a perturbation in mean position of the scattering centers 
and refractive index (n) in the focal volume of the focusing US transducer, which is the region 
of interest (ROI) for imaging. Coherent light picks-up a phase modulation from the ROI which 
is observed as an intensity modulation in the speckle formed through beating of exiting light. 
The modulation depth (M), obtained from the speckle intensity or field autocorrelation, is the 
measurement, from which both the optical (for example, absorption coefficient (\i a )) and me- 
chanical properties (for example, the shear modulus (G)) pertaining to the ROI are ascertained. 

The parent discipline whose extension UMOT is, is diffusing wave spectroscopy (DWS) 
wherein coherent light is employed to study multiple scattering media undergoing temperature- 
induced Brownian motion. The measurement is a correlation function (field- or intensity au- 
tocorrelation; g\(l) or g2(?) respectively) whose decay is used to quantitatively measure the 
mean complex modulus of elasticity of the medium. In UMOT we add a local, deterministic si- 
nusoidal perturbation to the stochastically driven multiple scattering medium at the focal region 
of the US transducer. This perturbation gives rise to a sinusoidal modulation on the correlation 
fluctuation. The speckle modulation mentioned earlier is a measure of the modulation in g\ (t) 
whose depth M is related to \l a and G of the US focal volume. 

A general assumption under which quantitative expressions connecting g\{l) to the proper- 
ties of the object is derived both in DWS and UMOT is that the optical wavelength A c is very 
small compared to the scattering mean-free-path (£ s ) [3, 4], which is the so-called weak scatte- 
ring condition. This would imply that correlation of light reaching the detector along different 
paths is negligibly small and only light traveling the same path contributes to Another 
assumption in regard to UMOT is that the phase increments (A0) introduced by US forcing of 
scattering centers within the ROI are uncorrelated. This is true when X a « £* where X a is the 
acoustic wavelength and t is the transport mean-free-path of light in the medium. In addition, 
one also assumes that the random process associated with the overall phase fluctuation picked- 
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up within the object is wide-sense stationary which is true when the object dimensions are large 
compared to X a . 

The above assumptions, particularly, the uncorrelatedness of A(j) within the ROI makes the 
phase measured in the g\ (t) modulation (i.e., the speckle interference term) a random variable 
with zero mean and therefore, carrying no useful information. On the other hand, when the 
acoustic excitation is of low frequency (say < 1 kHz) the assumption X a « I* is violated, and 
more importantly when the acoustic excitation is in the neighborhood of one of the dominant 
natural frequencies of vibration of the insonified US focal region, the movements introduced 
by the US forcing within a volume (£*) 3 are correlated and fall quite within one acoustic wave- 
length, and the corresponding phase increments picked up by light in that volume do not add 
up to a zero in the mean. The mean, (A(j)) av obtained from (along with M) enriches the 
UMOT data enabling recovery of additional parameters such as refractive index, n, and the 
scattering coefficient, \l s . For (A(j)) av to be nonzero, in addition to demanding that X a » £*, 
we should have a small enough ROI so that the intercept length of photon trajectory with the 
insonified ROI should have a good chance of being within £*. For this, our experience has 
been that if the ROI (which is usually an elongated hyperboloid with a narrow waist region) 
is bounded by 4 — 5 times (£*) 3 , (A(j)) av will have a consistent nonzero value. In addition, we 
have shown in Sec. 2 that the length scale involved in phase decorrelation is t n = t s ( 1 — (cos | ) ) 
which is larger than £* . 

Our aim in this work is to demonstrate by simulations and experiments that {A(j)) av measured 
from the speckle modulation can be nonzero if the above mentioned conditions are met. A 
brief summary of the rest of the paper is as follows: In Sec. 2 we put forth an expression 
for A(j) in terms of displacements of scattering centers and the local n. Here we also discuss 
in detail the conditions under which (A(j)) av is nonzero. Sec. 3 discusses the details of our 
simulation studies. Here we vary the volume of the ROI and X a (in effect the US frequency), 
and obtain histograms of (A0) measured at 256 x 256 elements in the CCD detector array. 
These histograms reveal a nonzero mean for (A0) when X a is large and the US frequency is in 
the neighborhood of a dominant natural frequency of vibration of the ROI, and the ROI volume 
is small. Sec. 4 describes an UMOT experiment wherein poly- Vinyl Alcohol (PVA) phantoms 
and chicken breast- tissue slabs are used as samples. From the histograms obtained (from the 
measured (A0) in all the elements in the detector array) we obtain (A(j)) av and demonstrate the 
sensitivity of the measured (A(j)) av to the local n in the ROI. Our concluding remarks are left to 
the last section (Sec. 5). 

2. Theoretical considerations 

Let N p be the number of scattering events experienced by photons in their passage through 
the insonified ROI. The scattering centers are themselves vibrating because of sinusoidal US 
forcing. The typical displacement of the j th scattering center about its mean position is given 
by 

Ar j(t) = u(r ; -)sin(k fl .r ; - - co a t). (1) 

Here u(ry) is the maximum displacement suffered by the j th scattering center, is the acoustic 
wave vector, rj is the position vector (at rest) of the j th scattering center and co a is the acoustic 
frequency. Fig. 1 shows a typical representation of scattering of light by a vibrating scattering 
center, with unit vectors in the incident and scattered photon directions (in a typical scattering 
event) represented by k{ and k s (k; and k? are the corresponding wave vectors) and the displace- 
ment by u (taken as along ki as shown in Fig. 1). The phase fluctuation experienced by the 
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photon following path p is given by [5] 



S$ p (t) = y £q j AY j (t). 



(2) 



Here q ; = k? — k ; . Since we consider electric field arriving along the same photon path alone 
to be correlated, the phase of the field autocorrelation of light reaching the detector along path 
p will have value: 



To find the average (A0(T)) flV which influences #i(t) we take the average of (A0 /7 (t))'s eval- 
uated over all the photon paths reaching the detector which have nonzero intercepts with the 
ROI, weighted by the probability density for the photon paths: 



Here s is the photon path length, and s\ and S2 are limits on path lengths s which pertain to 
those with nonzero intercepts with the ROI. The integral in Eq. 4 can also be evaluated, un- 




Fig. 1. Representation of a single scattering event from a point which suffered a maximum 
displacement of u. 

der weak scattering approximation, by first summing over (A$(t))'s pertaining to all paths 
passing through a particular scattering center and then summing the result over all the scat- 
tering centers in the ROI [5]. When this is done (A0(T)) flV can be nonzero only when there 
is scattering anisotropy for the medium to such an extent that (-T) 3 is large enough to have 
the intercept length of the photon trajectory with the ROI < I*. Moreover, while considering 
correlation of A0(t) one should use a scattering anisotropy factor g n = (cos|) to compute t 
instead of the usual g = (cos0). This results in a larger t (denote it by t n ) with which we can 
put limits on the volume of ROI so that (A0(T)) flV is nonzero (as demonstrated in Sec. 3). (This 
can be seen from the diagram of Fig. 1. The phase, A(j) picked-up in a typical scattering event 
is &o (k s — h) u : = kok.u = ko\k\ |u|cos [n — f ), where ko is the modulus of light propagation 
vector. Therefore the uncertainty in cos|, and not cos 6, determines the decorrelation of the 
overall phase fluctuation; and we use a modified g n = (cosf ) as the anisotropy factor to arrive 
at £,*). As we demonstrate through simulations, (A0(t)) <2V remains nonzero so long as the US 
focal volume V is bounded by 4 — 5 times (^) 3 . The reason for this persistence of correlation in 
(A0(T))av over a larger volume compared to (^*) 3 is that for the rather elongated shape of the 
ROI, the probability of a number of the intercept length of the photon trajectory with the ROI 
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(4) 



remaining within t n is large for ROI volumes considerably larger than (t n ) 3 . However, in our 
simulations we have used a cuboid ROI which is not elongated along any particular direction. 
The consequence of using this 'symmetric' object is that (A(j)(T)) av decays rapidly to zero as 
V increases from (i* n ) 3 to ~ 5(£*) 3 (See Fig. 2). In addition to the above there should also be 



■ 9 n = 0.97 



1 2 3 4 5 6 

Relative insonified volume with respect to (l* n ) 3 

Fig. 2. Variation of (A0(t)) (2V (arrived at through MC simulation) with the volume of the 
ROI. Beyond V = 4(£*) 3 , (A(j)(T)) av fluctuates rendering the average not reliable. 

'consonance' (i.e., directional bias along the US forcing) of the movement of scattering cen- 
ters in the insonified region introduced by the ultrasound forcing, for (A0(t)) gv to be nonzero. 
If X a « £*, there will be many 'periods' of US-induced oscillations within the insonified 
volume, which would mean that the resulting oscillations (or displacements) of the scattering 
centers will have equal propensity to be in a certain direction and its opposite. Therefore for 
'consonance' of the movement of scattering centers one should have X a » £*. This, how- 
ever, is not enough to ensure that the movements will have correlatedness in direction (and not 
randomly distributed) so that phase picked up is not averaged to zero. The jelly-like material, 
which is the model used to represent tissue, for the particular shape of the insonified volume, 
has natural frequencies of vibrations in the range of 100 — 400 Hz [6] corresponding to the 
axial translational mode. When the ultrasound forcing is at the natural frequency of vibration 
a resonance is observed in the oscillating movement of particles in the volume. In the vicinity 
of this resonant frequency, primarily only the prominent axial mode is excited (which sets up 
vibration along the US transducer axis) and the eigen values corresponding to the other modes 
such as non-axial translation, dilatation etc. are small. However when the forcing frequency 
is away from the resonant frequency of vibration, the input energy is more evenly distributed 
amongst all possible degrees of freedom with the result that oscillation amplitudes are more 
randomly oriented. Therefore when the forcing is at 1 — 2 kHz, even though one wavelength 
(~ m) generously covers the insonified region, owing to the apparent randomness of the appar- 
ent displacements the average phase picked-up is seen to be zero. 

Therefore, for the average phase to be nonzero we should have V ~ 4 — 5 times (t n ) 3 and 
h a » £% (for typical soft biological tissue, t ~ 1. 00 mm and t s ~ 0. 10 mm [7]), with the 
corresponding frequency close to the dominant natural frequency of vibration of the insonified 
region. In our experience (see Sec. 3) when the US frequency is below 1 kHz, i.e., X a ~ 1.5 m, 
(A(j)(T)) av evaluated over a cuboidal ROI of volume 3x7x3 mm 3 = 0.063 cm 3 is nonzero 
when g n ~ 0.97 and t n = 3.58 mm. If we take V = 5(t n ) 3 this works out to be 5 x 0.0461 cm 3 = 
0.2305 cm 3 , much larger compared to the volume of the ROI. 
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Using Eqs 3 and 4 one can obtain computable expressions for (A0(t)) qv in which the op- 
tical properties of the ROI and its size enter as parameters. In Sec. 3 we demonstrate through 
Monte Carlo (MC) simulations how (A0(r)) flv varies when (i) -^W is varied and (ii) Av the US 

forcing frequency is varied. We also demonstrate variation of (A(j)(T)) av when n is varied in the 
neighborhood of the ROI. 
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Fig. 3. Variation of (A0(t)) (2V with the ultrasound beat frequency. Beyond 1 kHz the varia- 
tion in the mean is large and the measurement is inaccurate. The error in the experimental 
measurements is in the range of 4 — 6%. 



3. Numerical simulations 

For numerical simulations we consider a cylindrical object of diameter 3 cm and height 5 cm 
(to represent PVA phantom), and a slab 3x3x5 cm 3 (to represent chicken breast sample). 
For PVA phantom, the optical properties used in the simulations are \i a = 2.85 x 10 -4 mm -1 , 
Us = 10.00 mm -1 , g = 0.89, and n = 1.35 [8], which is also varied to demonstrate the shift of 
(A(j)(T)) av . For chicken breast tissue (slab of dimension 3x3x5 cm 3 ) jl a , \l s and g are assumed 
to be 0.01 mm -1 , 23 mm~ l and 0.90 respectively in accordance with the properties of chicken 
breast tissue measured and reported in [9]. We use two identical focusing US transducers with 
focal regions intercepting in the object giving the ROI (for details, please see Section 4). The US 
focal volume of the individual transducer is assumed nearly hyperboloidal with axial length of 
~ 1.2 cm and diameter across the waist (the thinnest middle region) 0.25 cm. The intersection 
of two such focal volumes which is labeled the ROI is assumed to be a cuboid of volume V. 
In fixing these parameters we took recourse to computing the acoustic pressure (by solving 
the Westervelt equation as is done in [10]) inside the focal volume of the hemispherical PZT 
transducers that we, in fact, employed in the experiments. The US focal volume is fixed as the 
region where the pressure is greater than 100 th of its maximum pressure along the axis. The 
insonified volume, V, is varied in the simulations to demonstrate the dependence of (A0(r)) av 
on it. 

First we carried out two sets of simulations to demonstrate the effect of and Av on 

(A0(r))av m the first, through MC simulation, we compute the set of (A0(t)) <2V 's from which 
histograms of (A<j)(T)) av are plotted for V going from ~ 1.5(^*) 3 to ~ 6(^*) 3 . The results plotted 
in Fig. 2 show that the mean obtained from the histogram steadily goes to zero as V is increased 
from^l.5(^) 3 to^6(0 3 . 
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Fig. 4. Histogram of the phase fluctuation obtained from experiments and simulations. The 
object used is the PVA phantom. When the US frequency is (a) 1 MHz the mean is close 
to zero; when it is (b) 100 Hz the histogram has a clear bias with a peak around 140° (for 
experiment) and 145° (for simulations). The percentage error in phase extracted through 
measurements is 1 — 6% depending on intensity level. 



In the second set of simulations we varied Av from 100 Hz to 2 kHz in the steps of 250 Hz 
and computed (A^(t)) <2V . For Av > 1 kHz the average phase is close to zero, verified also for 
a high frequency of 1 MHz. Here, the dimensions of the ROI is kept at 3 x 7 x 3 mm 3 and the 
other optical properties are kept at values same as those used in the first set. The plot of the 
(A(f>(z)) av .vs. Av obtained through simulations is shown in Fig. 3, which also shows results 
from the experiments described later in Sec. 4. 



0.05r 0.12 




Fig. 5. Histogram of the measured phase fluctuation obtained from experiments and simu- 
lations. The object used is a slab of chicken breast which contains water in the ROI. When 
the US frequency is 1 MHz the phase fluctuation is random with the mean close to zero 
(a) and when it is 100 Hz the histogram has a clear bias with a peak around 110° (for ex- 
periment) and 105° (for simulations) (b). The percentage error in phase extracted through 
measurements is 1 — 8%. 

For arriving at (A0(r)) av we proceeded as follows: 5 x 10 7 photons are propagated using 
Monte Carlo simulation to compute the photon path distribution in the object in the transmis- 
sion geometry. We kept track of, primarily, the optical path length of photons in the insonified 
ROI and the photon weight from which the distribution of photon phase fluctuation is found, 
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using Eqs 3 and 4 at various exit points corresponding to the 256 x 256 elements in the CCD 
camera. (The details of MC simulation used are given in [11, 12]). From this we have plotted 
histograms of (A0(T)) flV for two US beat frequencies of Av = 1 MHz and 100 Hz. These are 
shown in Fig. 4(a) and (b) respectively for the cylindrical PVA phantom described earlier with n 
around the ROI assumed to be 1.33. Histograms obtained from experimentally generated phase 
distribution (details in Sec. 4) are also shown in the figure for comparison. When we repeated 
the simulations for the slab of chicken breast (with water, n = 1.33) in the ROI the histograms 
shown in Fig. 5(a) and 5(b) are obtained for Av = 1 MHz and 100 Hz respectively. It is seen 
that when frequency is 100 Hz the evaluated average phase in the two cases are 145° (for PVA) 
and 105° (chicken breast). When the simulation is repeated for chicken breast slab with glyc- 
erol in the ROI (n = 1.47) with frequency of 100 Hz we obtained the histogram shown in Fig. 
6 with (A0(T))av extracted as 120°. It is seen that the average phase fluctuation shifts when 
the n in the ROI changes. As stated earlier the above set of figures come with corresponding 
experimental results overlaid which will be explained in the next section. 




Fig. 6. Same as 5(b) but with glycerol injected around the ROI. It is seen that owing to the 
higher n of glycerol the average phase has increased to « 130° (for experiment) and 120° 
(for simulations). The percentage error in phase extracted through measurements is 1 — 7%. 



4. Experiments 

Our experimental set-up is similar to the one introduced in [13, 14] for UMOT which uses a 
source-locking technique to measure speckle modulation (of frequency as high as 1 MHz) with 
a slow 50 — 60 Hz CCD camera. For more details see [1 1] . It adapts the 4- or 3- intensity method 
of phase shift interferometry to compute the modulation and phase of interference patterns by 
recording 4 frames of intensities with preset phase differences between reference and object 
beams [15]. The system essentially uses a modulated laser beam locked to the driving oscillator 
of the US transducer. A low frequency acoustic excitation in the range of 100 Hz to 1 kHz is re- 
quired to obtain the resonance of vibrating ROI of jelly-like soft tissue-mimicking phantom [6] 
while US frequency of the order of MHz is demanded to achieve penetration depth of the order 
of em's in soft tissue. Here, we employed two identical, focusing US transducers (Panametrics, 
V394-SU-F50MM-PTF) operated in continuous mode with a small and adjustable frequency 
difference; one at V\ — 1 MHz and the other at V2 = V\ + Av (with the possibility of Av to 
be fixed accurately in the range 100 — 1200 Hz). The transducers are aligned, with the angle 
of inclination 60° (at the depth 1.5 cm inside the tissue samples), to have an intersecting com- 
mon focal region and US induced low frequency acoustic wave at the beat frequency Av at the 
region of intersection [11]. The experimental set-up is schematically shown in Fig. 7. During 
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Fig. 7. Schematic diagram of the experimental set-up. The CCD camera is focused onto the 
optical window. For a fuller description of the parts shown see Ref. [11]. 



the experiments, the samples are immersed in water (as shown in Fig. 7) for proper acoustic 
impedance matching. For entry and exit of light there are optical windows provided in the con- 
tainer. ROI is the narrow region of intersection of the two focal regions as depicted in Fig. 8. 
The near infrared (NIR) laser beam (of diameter 2 mm) of wavelength 785 nm is used to inter- 
rogate the US insonified ROI sent along the direction bisecting the angle subtended by the two 
US beams and the transmitted US modulated optical signal is collected using a CCD camera 
(JAI, CV-M10 BX, 25 frames/see) (see Fig. 8). The laser beam is intensity modulated whose 
injection current is a biased square wave of frequency Av and duty cycle of 25% composed 
of 40 — mA off-state and 60 — mA on-state, with Av equal to the US beat frequency. The laser 
diode is phase-locked to the function generator driving the US transducers and can be phase- 
shifted with respect to this latter drive signal. For a particular value of Av Hz we recorded four 
sequential speckle intensities with the CCD camera (with an exposure time of 40 msec each) 
each corresponding to phase shifts between the laser diode and US transducer drive signals of 
0, |, k and 3 1 radians. From these four intensities A(j) is computed using the formula given in 
Eq. A5 of Ref. [1 1] for each of the pixels, resulting in an ensemble of A0's. From this ensemble 
we compute a histogram of the distribution of phase in the range (—180°, 180°). 

In the first set of experiments we have used a cylindrical PVA phantom of dimensions and 
optical properties as given in Sec. 3. Experiments are conducted at number of frequencies start- 
ing at 100 Hz going up to 1 MHz, and the A(j) histograms are computed. From the histogram 
(A(j)} av is computed and plotted against Av in Fig. 3. For comparison, a similar plot obtained 
from numerical simulation is also shown. It is observed that the spread in the (A(j)) av values, 
in both the simulated and experimental results gets larger as the frequency increases. Beyond 
around 750 H z, (A(j)) av becomes both smaller and spread around the mean (Fig. 3). The ob- 
tained average phase is erroneous and unreliable beyond, say, 1 kHz- 

From the above experiment, the histogram obtained for two typical frequencies, namely, 
1 MHz and 100 Hz are shown in Fig. 4(a) and Fig. 4(b) respectively. For 1 MHz, A(j) is uni- 
formly distributed around 0° (without any bias towards any angle) indicating that (A(j)) av = 0 
(Fig. 4(a)), and for Av = 100 Hz it is peaked around 140° (Fig. 4(b)). In Fig. 4, we overlaid the 
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Fig. 8. Schematic diagram showing a typical cross section of the experimental set-up. Two 
focusing ultrasonic waves from the transducers are intersecting each other inside the tissue 
subtending an angle of 60°. A narrow laser beam (of diameter 2 mm) is incident onto the 
sample along the direction bisecting the subtending angle of the US focusing waves and 
the CCD faces the diode laser as shown to collect the transmitted light. 

corresponding simulated histograms obtained from Sec. 3 for comparison; and the match be- 
tween the simulated and experimental histograms and the averages is seen to be very good. This 
clearly demonstrates that the phase fluctuation has a nonzero mean when X a becomes very large 
compared to t n and the corresponding frequency is close to a natural frequency of vibration of 
the ROI. 

Similar histograms obtained by repeating the above experiment in the chicken breast slabs 
(with water, n = 1.33, in the ROI) are shown in Fig. 5(a) and Fig. 5(b) for 1 MHz and 100 Hz 
respectively. The average phase extracted from the low frequency histogram is seen to be 1 10°. 

In the next experiment we use an almost dry chicken breast slab with glycerol (n = 1.47) 
injected (using a hyperdermic syringe) so that the ROI is awash with glycerol. The new his- 
tograms (experimental as well as simulated) are shown in Fig. 6. It is seen that the (A(j)) av here 
is shifted to 130° for glycerol with n— 1.47 when compared to the first experiment for which 
the sample which is supposed to have water in the background with n = 1.33 gave the corre- 
sponding value as 110° .This clearly brings forth the role (A(j)) av can play in detecting changes 
in n deep within biological tissue samples. 

5. Conclusions 

In this study we demonstrate through simulations and experiments that the phase fluctuations 
carried by light propagating through the US focal volume can have a nonzero mean. Working 
under weak scattering approximation, the first condition to be met for the above mean to be 
nonzero is that the material should have a large g such that the US focal volume is bounded by 
4 — 5 times (^*) 3 . Another important requirement is that the phase picked-up by light scattering 
events within the ROI should not be averaged out by the randomness of the oscillations induced 
by the US forcing. For this it is not enough that X a » t n but also the frequency corresponding 
to the X a employed should be in the vicinity of the dominant natural mode of vibration of the 
insonified ROI. The nonzero mean of the phase of the ultrasound modulated light is an extra 
measurement that can be used to recover refractive index change in a localized region in the 
tissue specified by the US focal position. If in addition to the modulation depth M, phase is also 
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available through an UMOT experiment, albeit an average one, this can be used to construct an 
artificial hologram. This, in turn, can be used to recover a conjugate wave to back-propagate to 
the ROI forming either a \i a — or shear modulus image of the ROI. These can be the subjects of 
future investigations. 
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